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ABSTRACT 

In this paper we apply different methods to examine the possibility that a small group of 24 
asteroids dynamically linked to main-belt comet P/2006 VWiag, recently discovered by the 
Pan-STARRS 1 survey telescope, shares a common physical origin. By applying the Hierar- 
chical Clustering and Backward Integration methods, we find strong evidence that 1 1 of these 
asteroids form a sub-group which likely originated in a recent collision event, and that this 
group includes P/2006 VWisg. The objects not found to be part of the 1 1 -member sub-group, 
which we designate as the P/2006 VW139 family, were either found to be dynamically un- 
stable, or these are likely interlopers which should be expected due to the close proximity of 
the Themis family. As we demonstrated, statistical significance of P/2006 VW139 family is 
> 99%. We determine the age of the family to be 7.5 ± 0.3 Myr, and estimate the diameter of 
the parent body to be ^ 11 km. Results show that the family is produced by an impact which 
can be best characterized as a transition from catastrophic to cratering regime. The dynamical 
environment of this family is studied as well, including the identification of the most influen- 
tial mean motion and secular resonances in the region. Our findings make P/2006 VW139 now 
the second main-belt comet to be dynamically associated with a young asteroid family, a fact 
with important implications for the origin and activation mechanism of such objects. 

Key words: comets: general, comets: individual: P/2006 VW139, minor planets, asteroids: 
general, asteroids: individual: 300163 



1 INTRODUCTION 

Main belt comets (MBCs) are objects that are dynamically indis- 
tinguishable from main belt asteroids, but which exhibit comet- 
like activity due to the sublimation of volatile ice i Hsieh & JewittI 
L2OO6; Bertini 201 1). To date, six such objects have been discovered 
fasieh et alj|2012bi) . MBCs are important because they represent a 
new reservoir of comets in the Solar system, and may give insights 
into the role of main-belt objects in the primordial delivery of water 
to the Earth as well as provide constraints on the composition of the 
protosolar disk. 

Dynamical simulat ions performed by 

iFemandez. Gallardo, & Bruninil (l2002h show that it is un- 
likely that present-day objects with main-belt orbits originated in 
the Kuiper Belt, assuming the current configuration of the major 
planets. This indi cates that MBCs are likely native to the main 
asteroid belt. Still. iLevison et al] ( |2009|) show that during the early 
violent dynamical evolutio n of giant planet orbi t s, as required by 
the so-called Nice mode l jTsiganis et al.l |2003 : iMorbidelli et al.l 
I2OO5I : I Gomes et al.ll2005!) . some icy trans-Neptunian bodies may 
have been implanted in the asteroid belt. It is expected that most of 
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these objects are located in the outer belt, but they might be found 
anywhere with semi-major axes larger than about 2.6 au. 

Numerical simulations performed to assess the dynamical sta- 
bility of MBCs suggest that the orbits of the currently known ob- 
jects are stabl e, indicating that they are likely native to their cur- 
rent locations ('Haghighipoui^'2009t Blsieh et al.ll2012allbl) . The only 
exception is P/2008 Rl (Garrad d), which appears to be dynami- 
cally unstable jjewitt. Yang . & HaghighipouJ 120090 . These simu- 
lations did not take into account non-gravitational effects (e.g.. the 
Yarkovsky effect), however, that may play an important role in the 
long-term stability of the MBCs. 

Wherever MBCs originate from, their activity, although still 
not reliably unders tood, is likely triggered by the impact-excavation 
of subsurface ice jHsieh. Jewitt. & Femande3l2004 |jewitlll2012h 
because completely exposed surface ice is unstable against sub- 
limation at their heliocentri c distances over Gy r time-scales. Ac- 
cording to Hsieh ( 2009) an d lCapria et ai] ( l2012l) . this hypothesis is 
in reasonable agreement with the present-day impact frequency in 
the main belt. 

Thermal modelling by ISchorghofed fcOOSh and 



IPrialnik & Rosenberg! ( l2009l) shows that buried ice on a main-belt 
comet can in fact survive over the life of the solar system. How- 
ever, while thermal devolatization may not preclude the existence 
of present-day ice (and therefore MBCs) in the main asteroid 
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belt, there is the additional problem of coUisional devolatization. 
Each time an impact triggers activity in a MBC by exposing 
a small amount of subsurface ice to direct solar heating, that 
particular area of the surface is effectively devolatilized once 
that ice has sublimated away. Using estimates of the active areas 
required to produce the activity observed for MBCs, Hsieh (2009) 
determined that the observed activity was most likely tri ggered 
by a pproximately metr e- sized impactors. Using work by Cheng 
( I2OO4) and .Bottke et'al] ( l2005b() . it was then estimated that impacts 
by objects of this size in the main belt should occur roughly 
every 10* yrs. Over Gyr time-scales, the cumulative effect of 
such impacts could be to devolatilize a significant portion of the 
surface of an ice-bearing asteroid. More deeply buried ice could 
persist, safe from both thermal and collisional depletion, but would 
also therefore be inaccessible by activity-triggering impacts. Such 
asteroids would then be just as unlikely to exhibit activity in the 
present day and be identified as MBCs as other non-ice-bearing 
asteroids. 

Because of these reasons, it is important to understand 
the dynamical environment of each MBC and any possible 
links to coUisionally-formed asteroid families. Asteroid fami- 
lies are believed to originate in the cata strophic disruption of 
large asteroids (e.g. IZappala et al.l I2OO2I) . and are a fascinat- 
ing and challenging subject in themselves. A few tens of fam- 
ilies have been discovered to date across the main b elt (e.g. 
IZa ppala et al.ll995l : lMothe-Diniz. Roig. & Carvanol l2005^. provid- 
ing insights into the collisional history of the main belt (e.g., 
iMarzari! Farinella. & Davislll999t iBottke et alll2005ah . disruption 
event s over a size range inacc essible to laboratory experim ents 
(e.g., iMichel. Benz. & Richard son 2003; lOurda et al.l l2007h, the 
mineralogical structure of their parent bodies (e.g., Cellino et al.l 
I20021) . and many other subjects. 

Previous analyses show that three of the six currently known 
MBCs (133P/Elst-PizaiTO, 176P/LINEAR, and P/2 006 VW139) are 
assoc iated with t he larg e and old Themis family jHsieh & JewittI 
l2004lHsieh et alj20121]t) . which is thought to have formed from the 
catastrophic disruption of a ~400 km parent body 2.5±1.0 Gyr ago 
jNesvomv et al.l2003l) . Additionally, a fourth MBC (238P/Read) is 
close to the Themis family in orbital element phase space, and may 
have once belonged to the family (Haghighipour 2009). More sig- 
nificantly, 133P has been found to additionally belong to the small 
and young Beagl e family, which is th ought to have formed less 
than 10 Myr ago jNesvomv et al ]|2008h . This finding is significant 
because, using the devolatization rate derived by Hsieh (2009), if 
one assumes 133P is a primordial member of the Themis family, 
placing its age at ~2.5 Gyr, less than 10 per cent of 133P's sur- 
face would be expected to be unaffected by impacts by the present 
day. In contrast, 99 per cent of 133P's surface would be expected 
to remain unimpacted if the object is assumed to only be 10 Myr 
old (i.e., the age of the Beagle family). This p ossible link to youn; 
asteroid families is supported by the finding bv lHsieh et al.l ( l2012 
that P/2010 R2 (La Sagra) belongs to a small cluster that may have 
a common collisional origin. 

Comet-like activity in main-belt asteroid (300163) 2006 
VW139 was discovered by the Pan-STAR RS 1 survey telesco pe on 
Haleakala in Hawaii on 201 1 Novem ber 5 1 Hsieh et alj201ll) . with 
preliminary analysis of its activity llHsieh et al.ll2012b[) showing 
that it is likely to be a genuine comet (i.e., exhibiting activity due 
to sublimation), rather than a disrupted asteroid (i.e., whose ac- 
tivity is due to a recent impact) such a s P/2010 A2 (LINEAR) 
(e.g.. Ijewitt et alJl201Cl: ISnodgrass et al.l f20I0) and ( 596) Sc heila 
jjewitt et alj|201 iHBodewits et alJboTlh . ( iHsieh et aLll2012bl) also 
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Figure 1. Number of asteroids associated with P/2006 VW139 as a function 
of cut-off distance (in velocity space), expressed in m s~^. Two critical 
cut-off distance values are (i) 64 m s~^ when the number of dynamically 
associated objects increases significantly, and (ii) 70 m s^^ when the group 
around P/2006 VW139 merges with the Themis family. A potential sub- 
family within the Themis family is reasonably well defined in between these 
two critical values. 



examined P/2006 VWi39's possible links to dynamical asteroid 
fainilies and found that it is dynamically associated with the Themis 
family, but that it may also belong to a smaller sub-group that could 
represent a new previously-unidentified young asteroid family. 

In this paper, we conduct a detailed investigation into whether 
objects dynamically linked with P/2006 VWisg in fact have a com- 
mon physical origin. First, we identify members of the group using 
the Hierarchical Clustering Met hod (HCM). Then we use the Back- 
ward Integration Method (BIMi lNesvomv et al. 2002a) in order to 
further test the reliability of the group and to estimate its age. Fi- 
nally, we analyse the characteristics of this newly discovered family 
in greater detail. 



2 LINKING P/2006 VW139 TO DYNAMICAL FAMILIES 

We start our investigation by identifying the objects that are dy- 
namically associated with P/2006 VW139. The identification of as- 
teroid families is usually done in proper orbital element space using 
proper semi-major axes (a^), proper eccentricities (ep) and sines of 

f iroper inclination (sin(^p)) . For this analysis, we applied the HCM 
Zappala et alJll99"ol 1994) to analytically-de termined proper or- 
bital elements /Milan i & Knezevid 1990l . l 19941) available at the Ast- 
Dys web pagfl where we obtained elements for 398,841 asteroids 
(as of January 2012). 

The results obtained by the HCM are the same as those ob- 
tained by iHsieh et alj j2012bh . although here we use an updated 
proper element catalogue. Fig.[T]shows the number of asteroids dy- 
namically associated with P/2006 VW139 as a function of the cut- 
off distance, dc (in velocity space), the standard metric in the HCM. 
An intriguing sub-grouping appears at dc > 63 m s~^ before then 
merging with the Themis family at dc = 70 m s~^ . A complete list 
of asteroids that belong to this group is given in Table [T| Assuming 
that a cut-off value of dc = 70 m s^^ used by Nesvomy (2010) to 
characterize the Themis family is still appropriate, the sub-group 
around P/2006 VW139 is formally a part of this family. However, 
the latter result should be interpreted with some caution because the 
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Figure 2. The distribution of tlie members of Themis family (small crosses) 
and the sub-group around P/20()6 VW139 (filled circles) in the (ap,ep) 
plane. The most notable two- and three-body mean motion resonances lo- 
cated close to P/2006 VW139 group are marked as vertical dashed areas. 
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Figure 4. The stalactite diagram for 6374 Themis family members. Only 
the groups with at least Ncrit = 19 members are shown. Note that P/2006 
VW139 group is visible only in the one step, at 65 m s^^ (its location is 
emphasised by the circle). 
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Figure 3. The same as in Fig. |2] but in the (ap,sin(ip)) plane. Note that 
group around P/2006 VW139 is located on the edge of the Themis family. 



group occupies the outskirts of the Themis family. Also, one should 
bear in mind that the cut-off value appropriate for characterizing an 
asteroid family tends to decrease as the number of known aster- 
oids increases because the density of background objects becomes 
higher. Consequently, the appropriate value of dc for t he Themis 
family may be slightly smaller today than at the time of iNesvomvl 
( l2O10i) work. This is supporte d by the fa ct that list of the Themis 
family members produced by INesvomvl ((2010) does not include 
the asteroids from P/2006 VW139 group, although most of these 
objects were known at that time. 

On the other hand, as can be seen in Figs.|2]and[3] the P/2006 
VWi3g group is separated from the main part of Themis family 
by several two- and three-bod)|3 mean motion resonances (MMRs) 
and actually may be a real part of the Themis family simply sepa- 
rated by the MMRs. 

To visualize the structure of the whole Themis family, we gen- 
erate a stalactite diagram for the population of 6374 family mem- 
bers (Fig. |4ll. This type of diag ram is extensively used in other 
family identification papers (e.g.. lZapDala et al.lll990lll994ll995l : 

^ Three-body mean motion resonances are those whose resonant angle is 
defined as a linear combination of the mean longitudes of two planets and 



an asteroid iNesvomv & Morbidellilll998[) . 



iNovakovic. Cellino. & Knezevid [2oTlh . and is a highly effective 
means of visually displaying the structure of a given region. Here, 
however, we have an opposite example: a group, which we show 
later is likely real, is visible only in the one step used in the dia- 
gram, and would not be considered a family based on standard cri- 
teria for family identification. This situation is a good illustration 
of the limitation imposed to this approach in the region occupied 
by large and dense asteroid families such as Themis. 

The density of background objects in this region is high, and 
it is exceptionally difficult to separate nearby families due to the 
phenomenon of chaining0 Because of that, we are unfortunately 
unable to draw conclusive result whether this group belongs to the 
Themis family or not. However, this problem does not affect our 
further investigation of the P/2006 VW139 group. 



3 THE P/2006 VW139 FAMILY 

3.1 Membership and age from backward integrations 

To obtain a preliminary sense of the dynamical stability of asteroids 
that belong to the group, we calculate thei r Lyapunov times (Tiyap ) 
following the methodology described in lMilani~& Nobilil ( Il992h . 
Our results (see Table [T) show that the orbits of 19 objects can be 
considered stable (Tiyap > 100 kyr). Two unstable objects (221660 
and 2002 VU137) are likely trapped inside the 9J:4A resonance, one 
(51798) probably interacts with the 20J:9A resonance, and another 
(32868) appears to be destabilized by the 1J-I-3S — lA three-body 
MMR. The asteroid 60320 is located near, but outside the 20J:9A 
resonance (see Figs.llOlandll It. Its instability may be caused by an 
unidentified MMR of high order, or it may interact weakly with the 

^ A well-known drawback of the HCM based on the single linkage rule 
is the so-called 'chaining' phenomenon: first concentrations naturally tend 
to incorporate nea rby groups, thus forming a 'chain' (Zappala et al. 19941 : 
lEveritt et al.l200ll) . 
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Table 1. Members of the P/2006 VWiag Group 



Object 








sin{jpi£ 


H± 




Pv ± Cp^^ 


Df_ 




Clustering^ 


15156 


3.05542 


0. 


1647 


0.0437 


13.4 


0, 


,1032 ± 0.0207 


8.7 


151.5 


Y 


32868 


3.06753 


0. 


1639 


0.0467 


14.1 




(0.07) 


7.5 


45.3 


U 


51798 


3.05414 


0. 


1668 


0.0451 


14.8 




(0.07) 


5.5 


12.0 


U 


60320 


3.05300 


0. 


1602 


0.0398 


14.8 




(0.07) 


5.5 


91.7 


u 


117310 


3.04833 


0. 


1651 


0.0443 


14.7 


0. 


0587 ± 0.0084 


6.2 


3125.0 


12% 


144767 


3.05948 


0. 


1635 


0.0464 


15.8 


0, 


,0770 ± 0.0375 


3.3 


2325.6 


Y 


22075 1 


3.04909 


0. 


1608 


0.0433 


16.6 




(0.07) 


2.5 


1250.0 


Y 


221660 


3.03339 


0. 


1645 


0.0374 


15.2 




(0.07) 


4.5 


5.4 


U 


221760 


3.04065 


0. 


1621 


0.0407 


15.8 




(0.07) 


3.5 


3703.7 


Y 


250710 


3.04587 


0. 


1612 


0.0378 


15.9 




(0.07) 


3.4 


408.2 


6% 


255251 


3.04918 


0. 


1654 


0.0463 


17.0 




(0.07) 


2.0 


7142.8 


Y 


255333 


3.05631 


0. 


1626 


0.0489 


15.9 




(0.07) 


3.3 


970.9 


0% 


257688 


3.05167 


0. 


1629 


0.0519 


15.4 




(0.07) 


4.2 


4166.7 


Y 


261697 


3.05010 


0. 


1637 


0.0405 


16.5 




(0.07) 


2.5 


699.3 


Y 


294724 


3.06029 


0. 


1625 


0.0444 


16.8 




(0.07) 


2.2 


1428.6 


Y 


300163 


3.05350 


0. 


1613 


0.0379 


16.2 




(0.07) 


2.9 


2083.3 


Y 


305976 


3.04891 


0. 


1653 


0.0430 


16.5 




(0.07) 


2.6 


625.0 


0% 


2002 VUi37 


3.02924 


0. 


1627 


0.0384 


16.9 




(0.07) 


2.1 


4.8 


u 


2005 THi94 


3.05136 


0. 


1611 


0.0409 


16.6 


0, 


,0411 ±0.0165 


3.2 


6666.7 


0% 


2005 YCi49 


3.05162 


0. 


1613 


0.0510 


17.5 




(0.07) 


1.6 


4347.8 


Y 


2007 DPto 


3.04943 


0. 


1620 


0.0505 


16.8 




(0.07) 


2.2 


3333.3 


16% 


2007 UCi 


3.04305 


0. 


1643 


0.0418 


16.6 




(0.07) 


2.4 


5555.6 


18% 


2008 FGi4 


3.02182 


0. 


,1626 


0.0364 


17.4 




(0.07) 


1.7 


5555.6 


Y 


2011 UB54 


3.03709 


0. 


1625 


0.0378 


16.6 




(0.07) 


2.4 


800.0 


19% 



Proper semi-major axis, in a.u. 
^ Proper eccentricity 

Sine of proper inclination 

Absolute y-band magnitude, fr om the AstDvs website 
^ Geometric V-band albedo, from lMasiero et all bOl ih if available; otherwise the mean value for the group is given in brackets 

f Diameter, in km 

9 Lyapunov time, in kyr 

Notes about nature of clustering of secular angles, with "Y" indicating that clustering of secular angles exist, percentage values 
indicating the chance of clustering derived using clones, and "U" indicating unstable objects 



20J:9A resonance. We also note though that Lyapunov time of this 
asteroid is significantly longer than in the case of other unstable 
objects. More will be said about dynamical stability in Section [33] 

The fact that 19 out of 24 objects associated to the group 
around P/2006 VW139 are stable makes it possible to apply the 
BIM to search for possible clustering of secular angles in the 
past. As in the an alysis of the P/La Sagra cluster presented by 
iHsieh et alj j2012ah . the purpose of the BIM analysis here is two- 
fold: if suspicious clusterings in the longitude of the ascending 
node (n) and the argument of perihelion (ui) can be found, this 
analysis would not only allow us to estimate the age of the group, 
but would also provide support for the hypothesis that these objects 
share a common physical origin. 

As mentioned above, due to the proximity of the group around 
P/2006 VW139 to the Themis family, a distinction between the 
members of these two groups is very difficult to discern. There- 
fore, the probability of linking random interlopers with the P/2006 
VW139 group is much higher than in the case of "typical" fami- 
lies, where the fraction of interlopers is usually below 10 per cent 
( lMigliorinietalJl995l) . In some cases though, the fraction of inter- 
lopers may be as hi gh as ^25 per cent, as in the case of t he Eucharis 
asteroid family CNovakovic, Cellino. & Knezevidl2oTlh . 

As the members of the group cannot be reliably identified a 
priori, application of the BIM is complicated because the presence 



of interlopers could mask possible clustering that would detected if 
only real members were considered. Given that, we take a slightly 
unusual approach, where instead of searching for clustering that 
involves all asteroids linked to the group, we search for cluster- 
ings using sub-samples of potential family members. This method 
allows us not only to estimate the age of the group, but also to 
characterize the family by distinguishing individual members from 
interlopers Q We name this technique the Selective Backward Inte- 
gration Method (SBIM). 

Specifically, we perform the following four steps: 

(i) Starting with the list of dynamically linked group mem- 
bers, as identified by the HCM using the largest cut-off distance 
(dc = 69 m s^^) for which the group remains separated from 
Themis family, unstable asteroids are removed, i.e. those with 
Tiyap < 100 kyr. This allows us to search for possible clustering 
among as many asteroids as possible. 

(ii) Orbits of stable bodies are numerically integrated using the 
public domain 0RBIT9 software jMilani cfe NobilJll988h embed- 
ded in the multi-purpose OrbFit package. The dynamical model is 
purely gravitational and includes four major planets, from Jupiter 



* A similar idea was proposed bv lMarcus et al.l ( l201lb to search for fami- 
lies among Kuiper belt objects. 
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to Neptune, as perturbing bodies. To account for the indirect ef- 
fect of the inner planets, their masses are added to the mass of 
the Sun and the barycentric coiTection is applied to the initial con- 
ditions ( Milani & Knezevic 1992). Orbits of all objects are prop- 
agated over 20 Myr in the past. A online digital filtering is ap- 
plied to the output of numerical integrations in order to remove 
short periodic perturbatio ns. In this way, mean orbital elements 
jKnezevic & Milan? 2000) are obtained, which can then used to 
search for possible clustering of the secular angles. 

(iii) For the 19 stable asteroids in the P/2006 VW139 group, 
we check whether a clustering within 40° in at least one of the 
secular angles can be found for any sub-group with at least 10 
objects. The choice of the limit of 40° is based on previous re- 
sults obtained in application of the BIM jNesvomv etal J 120023 : 
iTsiganis. Knezevic, & Varvoglij2007l : lNovakovid2010l) . while the 
minimum number of objects is chosen assuming a maximum inter- 
loper fraction of 50 per cent. 

(iv) For each stable object that does not cluster, a set of 
orbit and "yarko" clones is generated t o check whether th e or- 
bit uncertainty or the Yarkovsky effect (lBottkeetal.ll2006l) may 
be the reason why clustering is not observed. In particular, 
ten orbit clones are drawn from 3-(t interval of the uncertain- 
ties of orbital elements listed in Table |2l assuming Gaussian 
distributions. Then, for each orbit clone, 10 yarko clones uni- 
formly distributed over the interval ±{da/dt)max are generated, 
where {da/dt)max is the maximum value of the semi-major 
axis drift speed due t o the Yarkovsky force. Using a model of 
the Yarkovsky effect jVokrouhlickvl il998l; IVokrouhlick^ 1 19991 : 



IVokrouhlickv & Farinellal 19991)^ and assuming therm al parameters 
appropriate for C-type objects JBroz & Vokrouhlickv 2008), we es- 
timate the maximum possible drift speed in the proper semi-major 
axis due to the Yarkovsky thermal effect, for a body of 1 km in di- 
ameter, to be {da/dt)max = 4 x 10^* au/Myr. In this way, a total 
of 100 clones of each object are produced. The orbits of all clones 
are numerically integrated in the framework of the same model as 
nominal orbits, but this time with the Yarkovsky force included in 
the modellfl Finally, a possible clustering of the secular angles of 
each of the clones, and that of the asteroids for which clustering of 
secular angles is observed in step (iii) is checked. The probability 
that the stable asteroids whose nominal orbits do not cluster should 
have exhibited clustering is estimated using the fraction of clones 
that exhibit clustering. 

Applying the above procedure, among 19 asteroids with 
Tlyap > 100 kyr, we find a sub-group of 11 objects that exhibit 
clustering in both and co about 7.5 Myr in the past. The clus- 
tering of Q is the only prominent one within the last 20 Myr and 
is very deep (within 32°). As such, we consider it to be a highly 
reliable event. 

On the other hand, although there are several clusterings of cj, 
none is particularly deep. This, however, is unsurprising, since we 
expect to find tighter clusterings of than of uj. The reason is that, 
in the region we consider, the secular frequenc}jf| gradient ds / dap 
is more than a factor of two smaller than d{g — s) /dap. Thus, small 
changes in a (e.g., due to the Yarkovsky effect) cause much quicker 
dispersions of lj in a way that is difficult to follow. Similar situa- 
tions were seen in the cases of the Veritas, Beagle and Theobalda 

^ For simplicity, the Yarkovsky effect is approximated in terms of a pure 
along-track acceleration, producing the same average semi-major axis drift 
speed da / dt as expected from theory. 

^ The s, g and g — s are frequencies of longitude of the ascending node ( fl), 
longitude of perihelion (zu) and argument of perihelion (lu) respectively. 




Time [Myr] 




Figure 5. The average of differences in the mean longitudes of the ascend- 
ing nodes f2 (top), and arguments of perihelion lu (bottom), for 1 1 dynam- 
ically stable asteroids that belong to the P/2006 VW139 group. The results 
are obtained in a purely gravitational model. The most important feature 
(visible in both angles) is clustering at 7.5 it 0.3 Myr ago, the likely age of 
the sub-group. 



families jNesvomv et al.l2002all2008l ; lNovakovi3l2010l) , for which 
the use of uj to search for clustering was ultimately impossible. 

Due to the aforementioned reasons, the remarkable match 
found between the deepest clusterings in SI and a; is a strong indica- 
tion that they are not the result of random fluctuations, but rather an 
indication that the objects considered share a common origin. We 
conclude that this group of 11 asteroids including P/2006 VW139 
is a real asteroid family formed 7.5 ± 0.3 Myr agoQ 

The last step in our analysis of clusterings of secular angles is 
to check whether orbit uncertainty or the Yarkovsky effect could be 
the reason why several stable asteroids linked to the group by the 
HCM do not cluster. This is done following the procedure described 
above in step (iv). The results of this analysis, presented in the last 
column of Table[T] show that this is unlikely for any of the objects, 
with probabilities below 20 per cent even in the most promising 
case of asteroid 201 1 UB54. Thus, only 1 1 asteroids can be reliably 
associated with the P/2006 VWiag family at this time. 



3.2 Statistical significance of the new family 

An important step in our analysis is to estimate a level of statistical 
significance of the P/2006 VW139 family. Because our conclusion 
that family is indeed real is mostly based on our findings of the 
clusterings of secular angles, we tested an hypothesis that these 
clusterings, simultaneously in both angles, can occur by a chance. 

In particular the following test was performed. In the same 
region as occupied by the P/2006 VW139 group, in terms of oscu- 
lating orbital elements (3.025 < a < 3.075 au, 0.12 < e < 0.21 
and 1.0 < i < 4.5°), we generated 2500 uniformly distributed test 
particles. For these particles we followed more or less the same ap- 
proach as for the members of the P/2006 VW139 group. Their orbits 



^ We use the phrase "P/2006 VWi3g group" to refer to the group of 24 
asteroids dynamically linked by the HCM, while "the P/2006 VW139 fam- 
ily" is used to refer only to those objects whose secular angles (f2 and ui) 
converge towards the same value. 
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Table 2. Stable members of the P/2006 VWisg group whose secular angles (H and a;) do not cluster 



Object 


oL 










mJ_ 
















117310 


3 0520756061 


1954589262 


3 67207 1 88 


89 08283896 

y .t./Oi.o JO y\j 


297 03599694 


1 62 03649209 




\J .V)\)\J\J\)\)\J<3\j i 


\.f .\j\j\.>\}\j\j 1 y^yj 


00001 ^{^f^ 


00094.470 


00094X^10 


0000^ 9 Q'^ 


250710 


^ ^4^^^^'^7S7l 

J . U'-tUOU J 1 / I 


U. 1 7 / Z JUU ill 


1 lQ^Af>(Mf^ 
i . / y j'^xjKfZ.yj 




1 S^^ Qf^l 70S04 


QO 79140Q9Q 
y\). 1 Lj'-\\jy I.y 




0.0000000937 


0.0000001605 


0.00001632 


0.00032530 


0.00032990 


0.00004392 


255333 


3.0533298729 


0.1298456130 


1.53037270 


298.13154200 


273.37048752 


284.24814758 




0.0000001916 


0.0000002877 


0.00001458 


0.00060860 


0.00062280 


0.00013090 


305976 


3.0625856084 


0.1226730090 


1.61519083 


215.05214592 


343.98604328 


221.68737391 




0.0000002950 


0.0000001935 


0.00001956 


0.00054650 


0.00055560 


0.00009934 


2005 THi94 


3.0455153624 


0.2008262940 


1.99615696 


191.88575916 


172.34882669 


110.98328885 




0.0000010430 


0.0000043270 


0.00003527 


0.00067090 


0.00288200 


0.00169800 


2007 DP70 


3.0516613516 


0.1343786668 


3.66093739 


145.64448060 


346.38036314 


4.96950838 




0.0000003648 


0.0000005328 


0.00004180 


0.00043090 


0.00055770 


0.00028460 


2007 UCi 


3.0382721562 


0.1874885082 


2.62463578 


165.94090018 


257.24768301 


275.35880634 




0.0000004401 


0.0000008379 


0.00002666 


0.00036510 


0.00062120 


0.00033210 


2011 UB54 


3.0330485525 


0.2025324904 


2.51664219 


36.42917884 


332.62952742 


52.94143685 




0.0000000996 


0.0000012870 


0.00005707 


0.00059980 


0.00071090 


0.00022170 



°' Osculating semi-major axis and formal uncertainty, in au 

^ Osculating eccentricity and formal uncertainty 
Osculating inclination and formal uncertainty, in deg 
Osculating longitude of ascending node, in deg 

^ Osculating argument of perihelion, in deg 

^ Osculating mean anomaly, in deg 



were integrated for 10 Myr and subsequently the Lyapunov times 
were determined. Then, 1900 test particles on stable orbits were 
selected and divided in 100 sets (each containing 19 particles). For 
these sets we checked whether clusterings within 32° and 55°, in 
n and w respectively, can be found for any smaller subset with 1 1 
objects. The characteristics of requested clusterings as well as the 
number of objects are the same ones as in the case of the P/2006 
VW139 family. Also, the fact that test particles are located in the 
same region of the main belt as the real family members, means 
that the dynamical conditions are the same, e.g. differential preces- 
sions of their secular angles are very similar. 

Performing these tests we did not find any appropriate clus- 
tering among 100 groups formed by the test particles. A clustering 
in o), within 55°, occurs at least once per 10 Myr in about 50 per 
cent of the cases. However, a clustering in Q., within 32°, is ob- 
served only in a few cases, and none of these were at same time as 
the clustering in lo. From this line of evidence we concluded that 
probability to find simultaneous clusterings in f2 and for any sub- 
group of 11 asteroids by the chance is less than 1 per cent. Thus, 
the statistical significance of the P/2006 VW139 family is greater 
than 99 per cent. 

3.3 Physical and dynamical characteristics of P/2006 VW139 
family 

To ascertain the physical characteristics of the members of the 
P/2006 VW139 group, we search for available data about their spec- 
tral types and albedos. Unfortunately, little is known about these 24 
objects. In particular, none of their spectral types have been de- 
termined to date. Geometric V^-band albedos, -pv, have been de- 
termined for four ob jects (Table [T), all from WISE observations 
jMasiero et al.ll201lh . The average geometric albedo derived from 
these four values is = 0.07. With their low albedos and prox- 



imity to the Themis family, most of these objects are likely C-type 
asteroids. The asteroid 15156 has a somewhat larger albedo than 
the other three bodies, but considering the uncertainty of this value, 
its albedo is still within the range expected for C-type asteroids. 
Spectroscopic observations are encouraged however to definitively 
confirm the spectral types of the members of the P/2006 VW139 
family, and characterize the family's level of physical homogene- 
ity. 

Next, we examine the dynamical characteristics of the P/2006 
VW139 group and the surrounding area. The positions of the most 
important MMRs have been already shown and discussed else- 
where in this paper. To better understand the role of these res- 
onances, and to check whether any connection between the dy- 
namical characteristics of the surrounding area and the apparent 
structure of the group can be established, we perform the following 
test: in a box slightly larger (3.018 - 3.077 au, 0.11 - 0.22, and 
0.9 — 4.6° in osculating semi-major axis, eccentricity, and incli- 
nation, respectively) than the one occupied by asteroids from the 
group, we generate 2000 randomly distributed massless test par- 
ticles. The orbits of these particles are then numerically integrated 
for 10 Myr, after which we calculate their synthetic proper elements 
dKnezevic & Milanil200d) and Lyapunov Characteristic Exponents 
(LCEs). 

The results are shown in Figs.|6]and|7] where we show colour 
coded maps of the obtained LCEs as functions of particle position 
in proper element space. We note that the figures indicate that most 
of the region surrounding P/2006 VW139 in orbital element space is 
stable, apart from two strongly chaotic zones centred at about 3.029 
and 3.075 au, corresponding to the 9J:4A and 11J:5A MMRs, re- 
spectively. Between these two zones, there are also some smaller 
areas that are weakly unstable. These weak instabilities are caused 
by higher order MMRs, such as the ones shown in Figs. 1 101 and 
[TTI Somewhat larger values of the LCEs are also found in the im- 



P/2006 VW139; A Main-Belt Comet Born in an Asteroid Collision? 7 



O 



0.18 
0.17 
0.16 
0.15 
0.14 
0.13 
0.12 



I 



140 

120 

100 

80 

60 

40 

20 





3.02 3.03 3.04 3.05 3.06 3.07 
Proper semi-major axis [au] 

Figure 6. The dynamical structure of the region ai'ound P/2006 VW139 
family. The location of main-belt comet P/2006 VW139 is marked by a 
star symbol. The colour scale indicates Lyapunov Characteristic Exponents 
(LCEs) (multiplied by 10'') for 2,000 test particles. The LCEs are inversely 
proportional to Tiy^p- For example, the value of 100 shown here corre- 
sponds to Tiyap = 10 kyr. 
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Figure 7. The same as in Fig.|6] but in the (ap,sin{jp)) plane. 
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Figure 8. The objects dynamically associated with the P/2006 VW139 
group plotted in the plane of proper frequency g of the longitude of per- 
ihelion versus proper frequency s of longitude of the ascending node. As 
in Figs. 1101 and 1111 circles represent group members whose secular an- 
gles (f2 and u)) cluster, while squares represent unstable objects. The lo- 
cation of main belt comet P/2006 VW139 is marked with the star symbol. 
Sizes of symbols are proportional to the estimated diameters of the cor- 
responding objects. The inclined yellow line indicates the position of the 
zi = g — ge + s — se secular resonance. The frequencies shown her e are 
calculated according to the synthetic theory l IKnezevic & Milanill2000t) . 
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Time [Myr] 

Figure 9. Time evolution of the resonant critical angle cr = ro — rog + Q — 
Slg of the zi = g — ga + s ~ sq secular resonance for asteroid 2008 FG14. 
The critical angle is in an anti-aligned libration state around 180° over the 
complete time period of 10 Myr covered by the integrations. 



mediate vicinity of P/2006 VW139. These appear to be caused by 
the 20J:9A and IJ-8S-1-IA MMRs. However, these resonances are 
weak and effective only in very narrow ranges of semi-major axis. 
An illustrative example is P/2006 VW139 itself, which is located 
very close to the 20J:9A resonance, but still does not show any sign 
of chaotic instability. Thus, we conclude that P/2006 VW139 family 
members are largely stable and unlikely to be significantly dynam- 
ically evolved. This conclusion is also valid for the single family 
member located on the opposite side of the 9J:4A resonance from 
the other family members, asteroid 2008 FG14. 

To develop as complete a picture of the dynamical environ- 
ment as possible, we also search for the presence of secular res- 
onances (up to degree six) in the region. Close to the group, we 
find only a non-linear secular re sonance zi = g — ge + s — sa 
jMilani & Knezevid 1 1991 1 19941) . Its location is shown in Fig. [8] 
Most of the objects from the group are far enough to avoid inter- 



actions with this resonance. This seems to be also true for asteroid 
2002 VU137, the only group member located on the other side of 
the g—ge + s—sa resonance (Fig. [8). However, due to the chaoticity 
of its orbit caused by the 9J:4A resonance, it is difficult to estimate 
its secular frequencies, g and s. Therefore, we are unable to draw 
firm conclusions as to the level of influence of secular resonances 
on this object. Finally, one family member, asteroid 2008 FG14, 
appears to be inside the g — g^ + s — sa secular resonance. 

To investigate this possibility, we plot the corresponding reso- 
nant critical angle a — vj — taJ(,+Q—^Ig in Fig.[9] The behaviour of 
(7, which oscillates around 180°, unequivocally confirms that 2008 
FG14 is locked within the zi secular resonance. Because of that, its 
association to the P/2006 VW139 family is less reliable, meaning 
that it could be an interloper 

In terms of proper semi-major axis, 2008 FG14 is located 
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significantly outside the assumed equivelocity curves shown in 
Figs.[TO]and[TT] The usual explanation for this would be Yarkovsky 
induced drift. However, the P/2006 VW139 family seems to be 
young, and the role of thermal effects is not so obvious. The 
maximum possible drift speed in proper semi-major axis due to 
the Yarkovsky effect acting on a 1 km body of {da/dt)max ~ 
4 X 10~* au/Myr (see Section l3Tt translates into a total drift of 
3 X 10~^ au over a period of 7.5 Myr As Yarkovsky induced drift 
is inversely proportional to the diameter of the body, for 2008 FG14, 
it reduces to the maximum value of 1.8 x 10^'' au since the family 
forming event. This value is almost one order of magnitude smaller 
than necessary to explain the observed displacement in Up. The ob- 
vious conclusion is that the Yarkovsky effect is not responsible for 
this object's anomalous displacement in Up, and that 2008 FG14 
could be an interloper in the P/2006 VW139 family. 

Insights into collisional physics may be obtained by analysing 
the inferred ejectio n velocity field (EVF) of the fragments (e.g., 
ICellinoetalJll999t) . An isotropic EVF can be computed using the 
well known Gaussian equations. These equations define ellipses 
in the (ap,ep) and (Op,sin(ip)) planes, whose shapes and orienta- 
tions are controlled by two angles, the true anomaly (/) and ar- 
gument of perihelion (uj) at the epoch of break-up. To connect 
these two ellipses with real family members, it is necessary to 
estimate the position of the family's barycentre (i.e., its centre 
of mass). Then, by adjusting the values of / and oj, and experi- 
menting with different ejection velocity changes (Av), it is pos- 
sible to find values of each of these parameters that best corre- 
spond to a given family, permitting us to infer some of the prop- 
erties of the EVF. However, it is well known that asteroid fam- 
ilies evolve over time with respect to their original post-impact 
configurations. They spread and become more dispersed due to 
chaoti c diffusion and gravitat i onal and non-gravita t ional perturba- 
tions (iNesvomv et al.ll2002bl: ICamiba et allboOSl ; iDeH'Oro et"an 
l2004l : lNovakovic. Tsiganis, & Knezevid20 1Cl). Young families like 
P/2006 VW139, however, may not have had enough time since their 
formation to be significantly affected by these effects, meaning 
that their ejection velocity fields (EVFs) may possibly be recon- 
structible. 

The first step in attempting to reconstruct the P/2006 VW139 
family's EVF is to estimate a position of the barycentre of the fam- 
ily. In order to do this, the absolute magnitudes (H; TableO of the 
members of the famil y are converted into diameters, D, using the 
standard relation (e.g. lBowell et al ]|l989l) 



0.175 



D (km) = 1329 



10" 



/Pv 



Measured albedos are used for the four objects for which these 
are available, while the average value of the measured albedos 
(Pv = 0.07) is assumed for all other bodies. Then, assuming the 
same density for all objects, we estimated their masses, which are 
then used to derive the position of the barycentre in proper orbital 
element space. Only the 1 1 objects classified as real family mem- 
bers are used in this calculation. The obtained barycentre position 



is at Op — 3.05480 au, Cp 



0.1644 and sin(ip) = 0.0431. The 



fact that asteroid 2008 FG14 may be an interloper has negligible 
influence on this result. 

As asteroid 15156 is the largest member of the family by far, 
the barycentre position is controlled by its position. As a conse- 
quence, the ellipses shown in FigsllOlandllllare shifted towards 
larger values of Up. This may be a result of a notably asym- 
metric ejection velocity field in terms of semi-major axis, caused 
by a large transverse velocity component. Such a situation of- 
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Figure 10. Members of the P/2006 VW139 group in the (ap,ep) plane. The 
superimposed ellipse marks an equivelocity curve, computed according to 
the equations of Gauss (e.g. Morbidelli et al. 1995), for the velocity change 
of Av = 60 m s~^, true anomaly / = 90° and argument of perihelion 
u] = 90° . Circular symbols represent group members for which Q and ui 
cluster about 7.5 Myr in the past, while squares represent unstable objects. 
Symbol sizes are proportional to the estimated diameters of the bodies plot- 
ted. Vertical shaded areas mark the approximate positions of local MMRs. 
The location of main-belt comet P/2006 VW139 is denoted by the star sym- 
bol. 



ten a r ises from cratering events ( Marzari et al.lll996t lNovakovi3 
|201(]I : IVokrouhlickv & Nesvomvll2oT H), though is sometimes ob- 
served in the cases of catastrophic fragmentation events as well 
( IZappalaetaLlll984l) . although these situations are less likely. To 
characterise the event which produced the P/2006 VW139 family, 
we estimate the size of the parent body by simply summing up the 
volumes of all 1 1 membero These volumes are obtained using the 
diameters listed in Table[T] Using this method, we calculate a nom- 
inal diameter for the parent body of Dps = 11.1 km. From this 
value we estimate the largest remnant to parent body mass ratio 
Mlr/Mpb to be 0.83. A reahstic value of Mlr/Mpb is likely 
somewhat smaller than the value mentioned here because the di- 
ameter of the par ent body estimated from known family members 
is a lower limit jPurda et~aL I l2007h . From this evidence, we can 
therefore conclude that the P/2006 VW139 family was not formed 
in a completely catastrophic fragmentation event. At the moment 
though, it is not possible to distinguish whether the formation of the 
family should be classified as a cratering or partially catastrophic 
event. In any case, the ratio Mlr/Mpb is high, and may be an 
explanation for observed asymmetry in the EVF. 

Alternatively, the calculated position of the barycentre may be 
wrong. It may be affected by exclusion of unidentified family mem- 
bers from our calculations. While the chance that any of eight sta- 
ble asteroids which do not exhibit clustering belong to the P/2006 
VW139 family is quite small, the status of the five unstable objects 
remains unclear Inclusion of these objects in the calculation of the 
barycentre position moves it to Op = 3.05604 au, ep = 0.1642 
and sin(ip) = 0.0437. In this case, the ejection velocity field be- 
comes even more asymmetric than before. The inclusion of these 
5 objects would additionally decrease the Mlr/Mpb ratio. This 
would make the explanation that the asymmetry of the EVF is a 
result of the impact event even less likely. Thus, we believe that at 



There are other, generally better, methods to estimate the size of a parent 
body ( Tanga et al. 1999; Durda et al. 2007.) , but we choose this approach 



here for simplicity. 
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Figure 11. The same as in Fig. [To]but in tlie (ap,sin(ip)) plane. The in- 
ner ellipse corresponds to an isotropic velocity field, while the outer one is 
produce assuming vertical component is 2.5 times larger than other two. 



ticularly around known MBCs, that may be previously unidentified 
young families. 

Meanwhile, now that the results of our application of the BIM 
indicate that the P/2006 VW139 group is likely a real asteroid fam- 
ily, it will be beneficial to provide observational confirmation of this 
dynamical result. Very little information about the spectral char- 
acteristics of asteroids belonging to this newly identified family is 
currently known. As such, we encourage physical characterizations 
of the members of this family, either photometrically via broad- 
band colours or spectroscopically, to help assess the plausibility of 
them being physically related in addition to be dynamically related. 
Furthermore, since family members are expected to share physical 
properties, the fact that P/2006 VW139 has been observed to exhibit 
cometary activity ap parently due to the sublimation of volatile ice 
( iHsieh et alj|2012bh . optical surveys for similar activity among its 
fellow family members may also be useful in the ongoing search 
for more MBCs. 



least some of the unstable objects associated to the group by the 
HCM are not real family members, but we are unable to draw a 
firm conclusion at this time. 

Finally, despite being currently fairly stable, asteroid 15156 
could have experienced some small displacement due to chaotic 
diffusion as it is very close to the slightly chaotic area centred at 
3.055 au (Fig.[6ll. Thus, it may have been less stable in the past, but 
then evolved onto its current stable orbit due to chaotic diffusion 
and/or the Yarkovsky effect. Even an exceptionally small amount 
of Yarkovsky induced drift of the semi-major axis would be enough 
to drive this asteroid from an unstable to a stable region (or vice 
versa). The maximum value of (da) 15156 is about 3 x 10"'^ au for 
7.5 Myr. Hence, the position of the barycentre might not be fully 
preserved. 

The last characteristic of the EVF that we want to draw atten- 
tion to is a notably larger spread in the inclination than in the ec- 
centricity (roughly by about a factor of 2). This could not be related 
to the position of the barycentre. Of course, it could be easily ex- 
plained by a large vertical velocity component of the EVF (see the 
two ellipses shown in Fig . lllb. Such situations are reported earlie r 
(e.g. IZappala et al J 1 19841 : iTsiganis. Knezevic. & Varvoglisll2007h , 
but to better understand the possible reasons for this, identification 
of additional family members is necessary. Hence, we postpone this 
analysis for a future work. 



4 SUMMARY AND CONCLUSIONS 

If the cluster of asteroids surrounding P/2006 VW139 does in fact 
have a common physical origin, this object would join 133P as the 
second MBC to be associ ated with a yo ung collisional family, sup- 
porting the hypothesis bv lHsiehldloog l) that such associations may 
be necessary for present-day sublimation on main-belt objects to 
occur. The lack of confirmed young family associations for any 
other MBCs does not necessarily undermine this hypothesis, as it is 
possible that many young families remain unidentified due to insuf- 
ficient numbers of their members simply not yet being discovered, 
or having insufficiently well-established orbits to perform dynami- 
cal linkage analyses. As more asteroids are discovered and orbits of 
known asteroids are improved by current and next-generation all- 
sky surveys (e.g. Pan-STARRS), it will be important to continue 
to search for small clusters of dynamically related asteroids, par- 
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